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Abstract 

A  numerical  technique  is  described  for  the  solution  of 
impact  on  finite  rods  in  which  the  strain  rate  is  assumed  to 
be  a  function  of  stress  and  strain.  Examples  are  given  for  the 
strain  rate  function  taken  as  a  power  of  the  dynamic  overstress 
and  the  results  are  compared  with  available  analytic  solutions 
and  approximate  uniform  strain  theory'. 
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1.  Introduction 

The  behavior  of  compression  impact  specimens  of  material  exhibiting  rate 
dependent  plastic  characteristics  is  of  interest  in  the  design  a,nd  interpretation 
of  experiments  for  investigating  dynamic  plastic  properties.  Previous  studies 
[1],  [2]  have  presented  approaches  by  way  of  analytical  solutions  of  the  partial 
differential  equations  describing  the  dynamic  deformation.  Analytical  solutions 
are  possible  only  if  the  general  partial  differential  equations  are  linearized, 
or  if  approximations  are  introduced  so  that  the  problem  is  reduced  to  solving 
ordinary  differential  equations.  Since  the  actual  behavior  of  structural  metals 
is  strongly  non-linear,  it  is  necessary  for  a  realistic  treatment  to  include  non¬ 
linear  relations  between  stress,  strain,  and  strain  rate,  as  well  as  irreversi¬ 
bility,  in  representing  the  plastic  behavior.  The  present  paper  describes  a 
scheme  of  a  numerical  solution  suitable  for  general  viscoplastic  laws.  The  tech¬ 
nique  is  illustrated  by  computations  for  two  types  of  viscoplastic  law.  More 
extensive  calculations  and  interpretations  with  experiments  will  be  given  in  a 
companion  paper. 

We  consider  the  problem  of  impact  on  a  finite  rod  as  shown  in  Fig.  1. 

A  rigid  body  of  mass  G  with  initial  velocity  parallel  to  the  axis  of  the  rod 
strikes  at  the  free  end  x  =  0  of  the  rod.  After  contact,  the  striking  mass  is 
assumed  to  stick  to  the  end  of  the  bar.  It  is  assumed  that  the  input  kinetic 
energy  GVo2/2  is  much  larger  than  the  maximum  elastic  energy  oq2AL/2E  which  the 
rod  can  absorb.  (Here  o and  E  are  respectively  the  static  yield  stress  and  the 
elastic  modulus  of  the  rod;  A  and  L  are  respectively  the  initial  cross  sectional 
area  and  the  length  of  the  rod.)  Hence,  the  elastic  strains  are  assumed  negli¬ 
gible  in  comparison  to  plastic  strains.  The  material  is  regarded  as  rigid-visco¬ 
plastic,  meaning  that  no  deformation  occurs  if  the  dynamic  stress  is  smaller  than 
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the  static  yield  stress  and  the  rate  of  deformation  is  governed  by  a  functional 
relation  between  strain  rate,  stress,  and  strain  if  the  dynamic  stress  exceeds  the 
static  yield  stress. 

We  introduce  the  following  dimensionless  notation: 


Table  1 


X 

O 

o 

4 

V 

a 

x  -  T 

L  - 

V  '  DL 

It 

pDL2 

o 

k  =  JL 

pAL 

o  e 
o 

Vo 

E1pD2L2 

\  - 

pD2L2 

Vo  "  DL 

°o2 

where 

X  =  distance  along  rod  from  impact  end 
T  =  time 

L  =  length  of  rod 

p  =  mass  density  of  the  rod  material 

G  =  mass  of  striking  body  (assumed  rigid) 

V  =  particle  velocity  at  X,T 

=  initial  velocity  of  striking  mass 

o  =  nominal  compressive  stress 

=  static  yield  stress  in  compression  (nominal) 

e  =  nominal  comppessive  strain 

D  =  constant  expressing  viscoplastic  behavior 

Ej  =  slope  of  strain  hardening  curve  (tangent  modulus  in 

plastic  range) 
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The  equations  of  dynamics  and  of  continuity  of  velocity  are  as  follows  in  the 
dimensionless  notation: 


s  =  -  v. 
x  t 


n.  =  -  v, 

t  X 


(1) 


where  a  subscript  denotes  partial  differentiation.  The  plastic  strain  rate 
behavior  will  be  expressed  by  (neans  of  a  function  F(s,n)  such  that 


nt  =  F(s,n)  >0  if  s  >  g(n) 


(2a) 


=  0 


if  | s |  £  g(n) 


where  s  =  g(n)  is  the  condition  for  plastic  flow  under  quasi-static  loading, 
with  g(0)  =  1.  F(s,n)  is  assumed  to  obey  the  following  conditions: 

F  >  0 
s 

F  <  0  (2b) 

n 

F(g(n),  n)  =  o. 


The  form  of  F  which  we  shall  consider  in  this  analysis  is 

n  =  (s  -  l  -  An)p  if  s  >  1  +  An 
t 

(3) 

=0  if  |s |  <  1  +  An 

where  p  and  A  are  constants  of  the  material.  The  power  p  corresponds  to  the  non¬ 
linearity  of  the  viscoplastic  behavior,  having  values  in  the  neighborhood  of  5 
for  several  structural  metals  (2).  The  coefficient  A  is  proportional  to  the  slope 
of  the  static  stress-strain  curve  in  the  plastic  range,  this  being  taken  as  linear 
for  simplicity.  For  a  more  general  form  of  F,  a  material  without  strain  hardening 
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can  be  characterized  by  F  =  0. 

n 

The  initial  and  boundary  conditions  of  the  problem  are: 


For  0  x  <  1: 


v(x,0)  =  n(x,0)  =  0 
s(x,0)  =  1 


For  0  £  t: 


v(l,t)  =  0  ;  v(0,0)  =  Vo  (4c) 

s( 0 ,t )  =  -  kvt(0,t)  (4d) 

The  boundary  condition  at  x  =  0  of  Eq.  (4d)  is  used  as  long  as  nt(0,t)  >  0.  When 
nt(0,t)  <  0,  there  is  an  unloading  boundary  £( t)  (see  Fig.  2)  along  which  nt(£,t)  =  0 
Then  Eq.  (4d)  must  be  replaced  by 

s(Cjt)  =  -  (k  +  ?)vt(c,t).  (4e) 

The  numerical  solution  of  Eqs.  (1)  and  (2)  subjected  to  conditions  (4)  will  be 
discussed  in  turn  in  the  following  sections  for  materials  without  strain-hardening 
(A  =  0)  and  with  strain-hardening  (A  >  0).  The  condition  (4b)  was  discussed  in  [1]. 

2.  No  Strain-Hardening 

For  materials  that  exhibit  no  appreciable  strain-hardening,  the  strain  rate 
function  F  can  be  assumed  independent  of  n.  Hence,  instead  of  Eqs.  (2)  we  write. 


nt  =  G(s)  >  0 


if  s  >  1 


if  s  <  1 


where  G  >0  and  G(l)  =  0.  Hence,  in  the  region  where  plastic  deformation  is 
s 


occurring,  Eqs.  (1)  and  (5)  reduce  to 


s  -  G  (s)s.  = 
XX  s  t 


0  if  s  >  1. 


(6) 
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This  is  a  quasi-linear  parabolic  equation. 

It  was  shown  in  [3]  that  the  solution  of  (6)  with  the  initial  and  boundary 
conditions  of  (4)  cannot  have  an  unloading  initiated  in  the  interior  of  a  plastic 
region.  The  proof  depends  on  the  maximum-minimum  principles  governing  equations  of 
the  type  of  Eq.  (6)  (see  [4]). 

When  Eq.  (6)  is  replaced  by  a  finite  difference  scheme  for  numerical  calcu¬ 
lation,  one  has  to  pay  attention  not  only  to  the  convergence  and  stability  of  the 
finite  difference  equations,  but  also  to  the  maximum-minimum  properties  of  the 
equations.  Otherwise,  physical  requirements  of  the  solution  may  be  violated.  The 
maximum-minimum  principles  for  various  finite  difference  equations  which  approximate 
Eq.  (6)  were  studied  in  [5].  The  results  show  that,  except  for  the  backward  differ¬ 
ence  method,  there  is  a  restriction  on  the  mesh  ratio,  which  depends  on  G  (s),  in 

s 

order  for  the  principles  to  hold.  If  we  denote  the  mesh  ratio  At/ (Ax) 2  by  r,  the 
results  of  [5]  show  that  the  weak  maximum -minimum  principle  holds  if 


r  <  G  (s)/2  for  the  forward  method 
s 

r  <  G  (s)  for  the  Crank -Nicolson  method 

s 

any  r  for  the  backward  method. 

If  Gg(s)  is  very  small  for  certain  s,  the  mesh  ratio  must  be  chosen  inconveniently 
small  unless  one  uses  the  backward  method.  For  instance,  without  strain-hardening, 
Eq.  (3)  is  written  as 


nt  =  (s  -  l)1'  if  s  >  1 
=  0  if  | s |  <  1 . 


(7) 


Hence  Gg(s)  =  p(s  -  1)'  For  p  >  1,  G  approaches  zero  as  s  approaches  one. 

Therefore,  in  order  to  use  the  same  mesh  ratio  for  the  whole  calculation,  we  shall 
use  only  the  backward  difference  method. 
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We  shall  derive  the  finite  difference  equations  directly  from  Eqs.  (1),  and 

represent  the  rod  for  this  purpose  by  a  model  consisting  of  discrete  mass  particles 

(Fig.  3).  Let  us  divide  the  rod  into  n  segments  (hence  Ax  =  1/n)  and  replace  the 

mass  of  each  segment  by  a  mass  particle  at  the  center  of  the  interval.  We  denote 

s.  n.  .  the  stress  and  strain  respectively  at  the  grid  point  x  =  iAx,  t  =  jAt, 

1  >1  1 » 3  '  .  . 

and  we  let  v.  .  denote  the  velocity  at  x  =  (i  -  i)Ax,  t  =  jAt.  The  strain  rate  will 
i  ,3 

be  written  as  i  .  instead  of  using  a  subscript  t.  Now,  Eqs.  (1)  and  (5),  together 
1  ,3 

with  the  boundary  and  initial  conditions  (4c)-(4e)  can  be  written  (see  Fig.  3)  as 
follows,  with  r  =  At/(Ax)2: 

Difference  equations: 


(a)  Impulse-momentum,  j  =  1,2... 


(8a) 


-n^k 

i  =  0:  s  .  =  -  (v  .  -  v  .  ) 

o,]  r  0,3  0,3-1 


i  =  1,2,  ...n:  s.  .  -  s.  ,  .  = - (v.  .  -  v.  .  n) 

i,3  1-1,3  r  1,3  1,3-1 


(b)  Viscoplastic  behavior,  j  =  0,1,2... 


(8b) 


i  =  0,1,2. ..n:  n .  .  =  G(s.  . )  if  s .  .  >  1 
i,3  i,3  i,3 


if  s.  .  <  1 
1,3 


(c)  Velocity-strain  rate,  j  =  1,2... 


0: 

v,  .  -  v  .  =  - 

1  . 
n 

1,3  0,3 

2n  0,3 

1,2. . .n-1: 

v. , ,  .  -  v.  .  = 

1  . 

—  n.  . 

i+l,3  i,3 

n  1,3 

n: 

1  . 

v  .  =  7T—  n.  . 
n,3  2n  1,3 

(8c) 
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Initial  Conditions  (j  =  0) 


i  =  0: 


v  =  v  ;  f)  =  2nv  =  G(s  ) 

o,o  o  o,o  o  0,0 


i  =  1,2. . ,n:  v.  =0;  n .  =0 

1,0  1,0 


s.  =  1 

1,0 


It  should  be  noticed  that,  with  Eqs.  (8b)?  the  unloading  condition  is  taken  into 

account  and  Eq.  (4e)  is  automatically  satisfied. 

One  can  eliminate  v.  .  from  Eqs.  (8).  If  this  is  done  one  obtains  for  the 

i,3  ■ 

general  case,  with  j  =  1,2,...;  i  =  l,2,...n-l: 


s.  ,  .  -2s.  .  +  s.  ,  .  =  —  [n .  .  -  f|.  .  ,] 
1+1,3  i-l,3  *  i,3  i,3-l 


(10a) 


=  -  [G(s .  .)  -  G(s .  .  )],  (10b) 

r  l,]  i,3-l 

making  use  of  Eq.  (5),  where  G(s.  . )  =  0  if  |s.  . |  s  1.  Equation  (10b)  corresponds 

1,3  i,3 

to  the  backward  finite  difference  form  of  Eq.  (6).  The  existence  and  uniqueness  of 

the  solution  of  s.  .  for  given  s.  .of  Eqs.  (10)  were  shown  in  [5]. 

1,3  i,3 -1 

For  given  v.  .  .,  Eqs.  (8)  can  be  solved  by  a  method  of  successive  approxi- 

i,3  -1 

mations  for  v.  .,  s.  .  and  n.  ..  By  assuming  an  approximate  value  v^°?  for  v  ., 

i,3i,3  1,3  0 ,3  °,3 

the  first  equations  of  Eqs.  (8a),  (8b)  and  (8c)  give  s^°\  rt^°^,  and  v^°?  respec- 

0,3  0,3  i,3 

tively.  The  second  equations  of  Eqs.  (8a),  (8b)  and  (8c)  then  give  s!j0?,  ru°?, 

1,3  1 ,3 

and  vi°] .  We  continue  this  process  until  v*'0?,  s^°?,  and  r/°?  are  obtained.  If 

2,3  n,3  n  ,3  n,3 

the  starting  value  VQ°j  is  the  correct  one,  then  the  last  equation  of  Eqs.  (8c) 

will  be  found  to  be  satisfied.  In  other  words,  if  corresponding  to  the  assumed 

(o)  ,(0)  , 

v  .  we  compute  ip .  ,  where 

0 , 3  3 

/0>SVM1.M  Ql) 

3  n,3  2n  n,3 


and  find  that  is  zero,  we  have  obtained  the  solution  for  time  t  =  jAt.  Other¬ 
wise,  we  assume  another  approximation  an<l  repeat  the  process. 
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Let  <5u  denote  the  error,  i.e.  the  difference  between  a  function  u  and  its 

.  +  .  (o) 

approximation  u  , 

r  (o) 

ou  =  u  -  u 

Then,  from  Eqs.  (8),  one  can  show  (cf.  [5])  that 

if  6v  .  >  0, 

o,D 

6\p  2  6v  .  2  6v  ,  .  i  .  i  Jv  .  >  0  (12a) 

n,]  n-1,3  0,3 

6s  .  <  6s  ,  .  <  .  <  6s  .  <  b. 

n,3  n-1,3  0,3 

Similarly,  if  6vq  .  <  0, 

6ip  <  6v  .  s  5v  ,  .  <  .  <  6v  .  <  0  (12b) 

n,3  n-1,3  0,3 


6s  .  i  Ss  ,  .  k  .  i  6s  .  >  0. 

n,3  n-1,3  0,3 

In  other  words,  the  absolute  values  of  the  error  in  the  calculation  of  v.  .  and 

s.  .  increase  as  i  increases.  The  errors  in  v.  .  are  of  the  same  sign.  Similarly, 
i»3  i,3 

the  errors  in  s.  .,  which  have  opposite  signs  to  v.  are  of  the  same  sign.  With 
i,3  i,3 

the  properties  expressed  in  Eqs.  (12)  there  is  no  difficulty  in  the  process  of 
successive  approximation. 

It  should  be  noticed  that  the  successive  approximation  can  be  improved  by 
considering  other  properties,  besides  Eq.  (12),  of  the  solution.  For  instance,  the 
fact  that  the  unloading  boundary  must  start  from  x  =  0  and  terminate  at  x  =  1 
implies  that 

if  s,  .  >  1 ,  s .  .  >  1  for  i  >  k . 
k,3  i,3 

Also,  it  can  be  shown  (see  [3])  that  for  1  <  j  <  N, 

v.  .  >  0  for  i  =  0,1, . n. 

i,3 
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At  time  t^  =  NAt  all  velocities  and  strain  rates  are  zero,  and  the  deformation  stops. 
The  strains  r).  .  are  obtained  from 

n. =  n.  .  +  n.  .At, 

which  is  essentially  a  backward  difference  expression  for  _. . 

For  illustration,  the  strain  rate  function  of  Eq.  (7)  is  taken.  In  the 
calculations,  p  is  taken  from  1  to  5,  Ax  =  0.05  and  At  =  0.025.  To  show  the 
accuracy  of  the  numerical  calculation  for  this  choice  of  mesh  size,  we  compare  the 
results  for  k  =  1,  v  =5,  and  p  =  1  with  the  exact  solution  of  the  differential 
equation  (see  [1]),  as  shown  in  Fig.  4.  Note  that  the  initial  condition  of  the 
problem  for  the  differential  equation  is  s(0,0)  =  ”,  while  that  for  the  finite 
difference  system  Eqs.  (9)  in  this  example  is  sq  q  =  201.  Nevertheless,  the  agree¬ 
ment  for  the  stress  distribution  is  fairly  good  at  t  .=  0.1,  which  is  only  four 
steps  from  t  =  0.  For  t  >  0.4,  there  is  practically  no  difference  between  the  exact 
and  numerical  solutions  for  both  stress  and  strain  in  the  rod.  The  time  of  unloading 

'  i  *  ,  •  »  ... 

t  is  1.665  and  the  total  time  of  deformation  t,.  is  2.120  in  the  exact  solution.  In 
o  f  • 

the  numerical  solution,  they  are  1.680  and  2.128  respectively. 

The  same  mesh  size  is  used  in  the  calculations  for  p  =  1  to  5,  with  the  same 
impact  data  k  =  1,  v  =  5.  Some  of  the  results  are  shown  in  Figs.  5  to  9,  Further 
discussion  is  deferred  to  later  sections,  after  the  presentation  of  the  solution 


for  strain-hardening. 


3.  Solutions  for  Strain-Hardening  Materials 

For  materials  that  exhibit  strain-hardening,  the  more  general  strain  rate 
function  Eqs.  (2)  must  be  used.  Although  the  system  of  equations  cannot  now  be 
reduced  to  a  second  order  parabolic  equation  in  s  alone,  we  shall  still  use  the 
backward  finite  difference  scheme  to  approximate  the  differential  equation.  This 
is  done  so  that  the  solution  will  exhibit  various  properties  found  for  non-strain- 
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hardening  materials,  when  the  material  has  very  small  but  finite  strain-hardening. 

We  use  the  same  model  as  before,  Fig.  3.  All  of  Eqs.  (8)  are  still  valid 
except  for  Eqs.  (8b)  which  must  be  replaced  by 
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0 


if 

if 


s . 


j 


(8b>) 


where  n.  .  =  n.  .  ,  +  n.  .At.  The  numerical  solution  of  Eqs.  (8a),  (8b')  and  (8c) 

i,3  i»3_i  i»3 

by  successive  approximation  is  essentially  the  same  as  that  illustrated  in  the 

previous  section.  The  only  difference  is  in  the  calculation  of  n.  .by  Eqs.  (8b’) 

when  s.  .is  known,  (n.  •  n,  of  course,  is  known.)  Eqs.  (8b')  for  s.  .  >  g(n.  .) 

1,3  i,3-l  i,3  i,3 

can  be  written  as: 


F(s.  . , 
i,3 


n. 


i,j-l 


.+  n.  -At). 

^  9  J 


If  this  equation  can  be  solved  for  n.  .  explicitly,  then  the  calculation  is  straight - 

i,3  " 

forward.  Otherwise,  an  iteration  scheme  may  be  devised  for  the  calculation  of  n.  . 

i,3 

from  s . 

i,3 

The  inequalities  expressed  in  Eqs.  (12)  are  still  valid  for  Eqs.  (8a),  (8b') 

and  (8c),  in  view  of  the  conditions  Eqs.  (2b)  imposed  on  F.  Therefore,  there  is  no 

difficulty  in  the  solution  by  successive  approximations  of  this  new  system  of 

equations.  The  existence  and  uniqueness  of  the  solution  follow  automatically. 

In  [3],  it  was  pointed  out  that  the  solution  of  the  differential  equation 

of  motion  with  strain  rate  function  represented  by  Eqs.  (2)  satisfies  the  minimum 

principle;  the  stress  s  cannot  have  a  minimum  value  in  the  plastic  region.  To  see 

that  the  finite  difference  equations,  Eqs.  (8a),  (8b1)  and  (8c),  do  possess  the 

same  minimum  principle,  one  can  eliminate  v.  .  from  Eqs.  (8a),  (8b')  and  (8c)  and 

1,3 

obtain  a  system  of  equations  similar  to  Eqs.  (10b): 
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s.  .  .  -  2s.  .  +  s.  .  .  =  —  [F(s.  .,  n.  •)  -  F(,s .  .  ,  n.  .  )] 

i+l,]  i,3  1-1,3  r  i,3  i,3  i,3-l  i,3-l 


(13a) 


where  i  =  1,2, . n-1.  The  equations  for  i  =  0  and  n  are  not  needed.  Now,  by 

the  mean  value  theorem,  the  right-hand  side  of  Eq.  (13a)  can  be  written  as 

1,3  1,3  i,3-l  1,3-1  i,3  1,3-1  s  i,3  i,3"l  n 

=  (s.  .  -  s.  .  . )  F  +  n.  .  F  At 

1.3  1,3-1  s  1,3  n 

=(s.  .  -  s .  .)F  +  F(s .  .,  n.  .)  F  At 

1.3  i»3“l  s  i,3  1,3  h 

where  F  =  F  (g.  .,n.  .),  F  -  F  (i.  .,  rj.  .) 

s  sv“i,3  i,]'*  n  n  i,3  i,3 

s.  .=0s.  .  +  (1  -  0 )s .  .  , 

i,3  i,3  1,3-1 


n.  .  =  0n.  .  +  (l  -  0)n.  .  . 

1,3  1,3  1,3-1 


0  <  0  <  1. 


Therefore , 


F(si,r  ni,j)  = 


and  Eq.  (13a)  becomes 


1-F  At 

n 


>  .  .  -I  ,  q*  .  *1^"^  n  s .  .  ^)  f  1]  , 

i,3  ~±  i  ,3  “I  i,3  1,3-1  s 


s .  .  .  -  2s.  .  +  s.  .  .  =  -  [AtF  F(s .  .  11...)  +  (s.  .  -  s.  .  .  )F  ]  • 

1+10  1.:  1-1.3  rd-F  it)  ”  1'3-1’  's 

n 

(13b) 

Since  Fg  >  0,  <  0  by  Eqs.  (2b),  Eq,  (13b)  satisfies  the  conditions  of  the  minimum 

principle  as  shown  in  Lemma  II  of  [5].  Thus  s.  .  cannot  have  its  minimum  at  an 

i,3 

interior  point. 

As  an  illustration,  we  take  the  strain  rate  function  represented  by  Eqs.  (3) 
with  X  =  1.  The  results  for  k  =  1,  v  =  5  and  p  =  1  to  5  are  shown  in  Figs.  6  to  10 . 
In  all  these  calculations,  the  same  mesh  size  (Ax  =  0.05,  At  =  0.025)  is  used  as  for 
the  previous  example. 
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For  p  =  1,  the  problem  is  linear  and  the  analytic  solution  can  be  obtained 
by  using  the  Laplace  transform.  This  is  shown  by  the  bold  face  lines  in  Fig.  7B 
for  the  solution  before  unloading  occurs.  It  is  seen  that  the  differences  between 
the  exact  and  numerical  solutions  are  indistinguishable  for  larger  t. 

To  investigate  the  accuracy  for  p  other  than  one,  we  consider  the  example 
for  p  =  3.  When  Ax  is  decreased,  the  solution  does  not  seem  to  improve  in  any  way. 
For  smaller  At,  however,  the  solution  does  depend  on  At  as  shown  in  Fig.  5.  Here 
we  see  that  the  unloading  time  tQ  and  the  total  strain  n(0,  t  )  at  the  free  end  vary 
nearly  linearly  with  At.  For  more  detailed  comparisons,  the  solution  for  p  =  3 
with  At  =  0.0125  (which  is  half  the  original  interval  At  =  0.025)  is  shown  in  Figs. 
6B  and  7B  by  the  dotted  lines.  The  results  are  fairly  good  for  x  t  0.  In  general, 
the  result  is  qualitatively  satisfactory  even  for  x  =  0,  if  not  quantitatively. 

The  difficulty  in  the  accuracy  of  the  solution  at  x  =  0  is  the  fact  that  the  ori¬ 
ginal  differential  equation  has  a  singular  initial  condition  while  the  approximate 
finite  difference  equation  does  not  (and  from  the  numerical  point  of  view  cannot) 
have  this  kind  of  singularity. 

Figures  6-12  present  curves  showing  plots  of  various  quantities  as  functions 
of  x  or  t.  Figures  6-10  show  results  for  impact  conditions  k=l,  v=5;  while 
Figs.  11  and  12  are  for  k  =  1  and  10,  with  vq  ranging  from  0.1  to  20. 

The  most  interesting  single  feature  of  these  results  is  perhaps  the  appear¬ 
ance  of  qualitative  as  well  as  quantitative  contrasts  between  the  cases  X  =  0 
and  X  =  1,  i.e.  between  a  pure  viscoplastic  material  (perfectly  plastic  behavior 
with  yield  stress  dependent  on  strain  rate),  and  a  modified  form  of  viscoplastic 
behavior  in  which  strain  hardening  is  considered.  For  example,  considering  final 
strains  plotted  as  functions  of  x,  Fig.  8A  (for  X  =  0)  shows  monotonically  de¬ 
creasing  functions  whereas  in  Fig.  8B  (for  X  =  1)  the  curves  for  p  >  z  show 
pronounced  minima;  that  for  p  =  1  monotonically  increases.  The  curves  in  Figs.  11 
and  12  (for  a  variety  of  k  and  vq  magnitudes)  also  exhibit  characteristic  features 
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of  being  concave  up  for  X  =  0  and  concave  down  for  X  =  1,  at  the  point  x  =  1. 

This  difference  between  the  two  media  was  proved  analytically  for  p  =  1  in  [3], 
and  the  present  calculations  suggest  that  it  holds  also  for  the  non-linear 
materials  considered  here. 

The  curves  of  Figs.  11  and  12  show  final  strains  n^(x)  in  ratio  to  the 

ii  t 

strains  computed  from  an  analysis  which  assumes  that  strains  and  stresses  are 
at  all  times  uniform  along  the  specimen  rod.  This  simplified  analysis  is  outlined 
in  the  next  section.  It  is  seen  that  the  final  strains  tend  toward  the  result  of 
the  uniform  strain  theory  as  either  k  or  vq  increases ,  but  the  effect  of  increasing 
k  is  particularly  strong. 

Further  remarks  on  physical  interpretation  of  these  results  are  given  in 
the  final  section. 

4.  Uniform  Strain  Analysis 

As  the  striking  mass  or  impact  velocity  is  increased,  the  strains  tend  toward 
a  uniform  distribution  along  the  rod.  An  approximate  theory  assuming  uniform  stress 
and  strain  is  presented  here,  appropriate  for  sufficiently  large  k  or  vq. 

If  v(t)  is  the  velocity  of  the  striking  mass,  since  all  quantities  are  inde¬ 
pendent  of  x  the  equations  of  motion  and  strain  rate  behavior  are 

s(t)  =  -  kvt(t) 

nt(t)  =  F(s,n)  =  v(t) 

Eliminating  s  and  v  between  these  equations, 

nt  =  F[-kntt,  n] 


(14a) 
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with  the  initial  conditions 


n(o)  =  o,  nt(o)  =  vq. 

Equations  (14)  can  be  reduced  to  a  first  order  equation.  Let 


(14b) 


Then 

and  Eqs.  (14)  become 


=  q. 


'tt 


= 

dn 


q 


q  =  F[-kq  n] 

q  =  v  when  n  =  0 . 
o 


(15a) 

(15b) 


If  Eq.  (15a)  can  be  solved  for  dq/dn,  it  is  convenient  to  put  it  in  the  form: 


—  =  f(q,n).  (16) 

Conventional  methods  of  integrating  first  order  differential  equations  can  be 
applied  to  Eq.  (16)  with  initial  condition  (15b).  For  comparison  with  the  results 
of  previous  sections,  we  take  F  as  expressed  by  Eq.  (3).  Then  Eq.  (16)  becomes 


dn 

dq 


kq 

1  +  Xn  +  q1^ 


(17a) 


with  the  initial  condition 


q  =  v  ,  n  =  0.  (17b) 

o 

Our  object  is  to  find  the  final  strain  n^  which  is  the  value  of  n  at  q  =  0i 

A 

The  integration  of  Eq.  (17)  was  done  by  electronic  computer,  and  the  values  of  nf  have 
been  used  in  presenting  the  results  of  the  complete  solution  in  Figs.  11  and  12. 
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Equation  (IV)  can  be  solved  analytically  in  a  closed  form  when  X  =  0 
with  p  integer  or  X  i  0  with  p  =  1.  For  X  =  0  with  p  integral,  the  solution  is 
given  in  [2].  For  X  /  0  and  p  =  1,  the  solution  is 


ft 

Xnf  + 


l 


( 1  +  ,  /?. u )  exp 


-  /2y 
1  +  /2 v 


if  v  2  =  8y 
o 


(18) 


if  v  2  <  8y 
o 

where  u  =  —  kv  2. 

2  o 

5.  Discussion 

In  the  preceding  work  the  non-dimensional  variables  listed  in  Table  1  were 
used  for  convenience  in  the  computations.  For  physical  interpretations,  these  have 
disadvantages  of  mixing  physical  quantities  specifying  the  impact  conditions  with 
those,  describing  the  material  behavior.  Also  the  quantities  specifying  plastic  and 
viscous  properties  are  intermingled.  For  example,  the  impact  data  are  specified  by 
the  quantities  k  =  G/pAL  and  vq  =  Vo/DL,  and  the  parameter  X  =  E1pD2L2/oo2  is  used 
to  specify  the  strain  hardening, behavior .  These  are  convenient  for  numerical  work, 
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but  make  it  difficult  to  draw  comparisons  between  different  materials  and  impact 
conditions.  The  following  dimensionless  parameters  have  more  obvious  physical 
meaning  and  are  more  convenient  for  physical  interpretations: 


In  the  new  parameters  of  Table  2,  the  quantity  /a  /p ,  with  dimensions  .of 
velocity,  is  a  convenient  reference  characteristic  of  the  material.  There  are 
two  independent  parameters  specifying  the  impact  data ,  which  can  be  defined  in 
various  ways  but  in  any  case  enable  one  to  specify  the  initial  energy  and  initial 
momentum  of  the  striking  mass.  The  most  useful  impact  parameter  is  e°,  which 
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represents  the  final  strain  in  the  specimen  under  the  idealized  conditions  of 
uniform  strain  and  stress,  all  the  energy  being  imagined  as  absorbed  at  the  yield 
stress  ao  with  neither  strain  hardening  nor  rate  sensitivity .(  Since  e°  is  also  the 
initial  kinetic  energy  of  the  striking  mass,  in  a  non-dimensional  form,  the  energy 
ratio  R  of  initial  energy  to  the  maximum  elastic  strain  energy  can  be  written  as 


R 


GV  2 
o 

2 


2E  o 


(19) 


A  necessary  condition  for  validity  of  a  rigid-plastic  treatment  is  that  R  >>  1; 
for  the  impact  problem  this  would  be  also  a  sufficient  condition.  For  e°  of  the 
order  of  0.05  or  more,  R  would  be  roughly  100  or  larger  for  one  of  the  structural 
metals,  and  presumably  a  rigid-plastic  treatment  would  apply  with  good  accuracy. 

Similarly,  the  parameter  t°  is  the  (non-dimensional)  duration  time  of  the 
deformation;  it  is  also  the  initial  momentum  in  a  non-dimensional  form. 

Estimates  of  the  approximate  ranges  of  numerical  values  of  /aQ/ p  and  other 
parameters  representing  material  properties  were  given  in  [2]  as  follows: 


o  (ksi) 
o 

/aQ/p  (in- sec  "b 

E  /a 

P  o 

=  6 

P 

D 

=  vL  ( in . ) 

Mild  steel 

30 

6,500 

6-9 

■ 

25-200 

Aluminum 

alloy 

6061-T6 

40 

12,500 

1-4 

4-8 

0.05-2 

With  these  ranges  in  mind,  the  physical  significance  of  the  illustrative 
cases  treated  in  this  report  can  be  seen.  Results  have  been  given  for  k  =  1,  10; 
v  =  1,  5,  20;  A  =  0,  1.  Effects  of  changes  of  one  of  the  impact  variables  k  or 
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e°  (or  4>q) ,  or  of  the  material  parameters  v  or  0,  can  be  observed  if  all  the  remaining 
quantities  are  held  constant.  For  example,  suppose  v  and  8  have  constant  values; 
changes  of  k  or  of  vq  then  show  the  effects  of  changing  impact  conditions  on  specimens 
of  a  given  material.  Table  3  shows  the  numerical  values  of  <j>o,  e°,  and  8  assuming 
v  =  20,  which  correspond  to  the  magnitudes  of  k,  vq,  and  X  used  in  the  calculations 
of  this  report. 


Table  3 

Interpretation  of  impact  data  with  v  =  20 
e  =  n/v2  =  n/400 
8=0  for  X  =  0 
6  =  400  for  X  =  1 


V 

4> 

£f 

o 

o 

il 

A! 

k  =  10  ■ 

0.1 

0.005 

.125xl0"4 

1.25xl0-4 

1 

0.05 

.125xl0~2 

.0125 

5 

0.25 

.0313 

0.313 

20 

1.0 

.5 

5 

The  results  shown  in  Figs.  11  and  12  thus  correspond  to  tests  on  a  material 
of  moderately  strong  rate  sensitivity  and  either  no  strain  hardening  or  very  large 
strain  hardening.  The  range  of  impact  energies  considered  is  too  wide  to  be  realistic 
at  the  small  values  of  vq  the  condition  R  >>  1  for  a  rigid-plastic  treatment  could 
hardly  be  satisfied,  while  at  the  largest  values  of  v^  and  k  the  strains  are  unrealis¬ 
tically  large. 
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Alternatively  one  can  take  the  impact  conditions  constant  and  consider  the 
effects  of  changing  material  parameters.  However  the  relation  3  =  Av2  shows  that  if 
A  is  constant,  a  change  in  3  accompanies  one  in  v.  The  present  calculations  are 
therefore  not  suitable  for  comparing  the  responses  of  different  materials  under 
constant  impact  conditions.  Further  calculation  for  this  purpose  has  been  carried 
out,  and  will  be  presented  in  a  companion  report  [6]. 
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